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Abstract 

The fluid models mentioned in the title are studied in a modified 
approach, based on two formulas for the mass function. All character- 
istics of the fluid are expressed through a master potential, satisfying 
an ordinary second order differential equation. Different constraints 
are imposed on this core of relations, finding new solutions and de- 
riving the classical results for perfect fluids as particular cases. All 
charged anisotropic solutions, all conformally flat and all uniform den- 
sity solutions are found. A large class of solutions with linear equation 
among the two pressures is derived, including the case of vanishing 
tangential pressure. The mechanism responsible for the appearance of 
equation of state is elucidated. 

1 Introduction 



The description of gravitational collapse and evolution of compact objects 
under various conditions remain among the important problems of general 
relativity. They are described by spherically symmetric relativistic fluid mod- 
els where the metric depends on time and radius. In general, they possess 
shear, expansion and acceleration, which makes them rather hard to study. 
In many physical situations, however, shear-free perfect fluid models are a 
good approximation. The history of their study is long and rich in rediscov- 
eries p] . McVittie used an ad hoc ansatz to find a solution [2] . The solution 
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of the Einstein equations may be reduced to a single ordinary second order 
differential equation [3], [1]. Its integration may be done by Lie point symme- 
try methods [1], [5], [6], [7J, [8]. Another approach is to search for solutions 
which have the Painleve property [9]. Other solutions have been found too 

The study of charged shear-free perfect fluids started with the papers of 
Shah and Vaidya [UJ and Faulkes [T2]. The master equation acquires an 
additional non-derivative term. The three approaches mentioned above were 
applied to the charged case too finding new solutions [13] , [H] , [15] , [16] , [17] , 
[T8] , [19] , [20] , [2TJ . It was noticed that both in the neutral and in the charged 
case the second Weyl invariant is tightly bound to the arbitrary functions in 
the master equation. Solutions which obey an equation of state (EOS) have 
been discussed in the neutral [3J, |22j and the charged case [23], [21], [25J, 
[26] and are known as the Wyman and the charged Wyman solution. 

All perfect fluid solutions with uniform energy density have been found 
by Kustaanheimo [lj, [27J. There are two branches - one in elementary 
functions and one given implicitly. Different properties of these solutions 
have been discussed [2S],i2n],[3^,[3TJ,[3J,[33],[r6],[31]. In some of these 
papers it was noticed that conformally flat solutions are particular cases 
of uniform density solutions. Charged shear-free perfect fluid spheres with 
uniform density are necessarily static [35] . As for the conformally flat perfect 
fluid solutions, all of them have been found by Stephani, using an embedding 
in a flat sixdimensional space [TJ, [36J. Some approximate conformally flat 
solutions for anisotropic fluid with slow motion were studied recently [37], as 
well as conformally flat solutions for radiating perfect fluids [38], [39]. Close 
to them are solutions admitting conformal symmetry [lj, [3Dj, [41] . 

On the other hand, different mechanisms have been identified through the 
years which create pressure anisotropy in stellar models and make the fluid 
imperfect [12]. Such are the exotic phase transitions during gravitational 
collapse, the existence of a solid core or the presence of a type P superfluid. 
Viscosity may also be a source of local anisotropy as well as the slow rotation 
of a fluid. It has been shown that the sum of two perfect fluids, two null fluids 
or a perfect and a null fluid may be represented by effective anisotropic fluid 
models [43J. Recently it was pointed out that the same is true for perfect 
fluids with charge, bulk and shear viscosity [H]. 

The aim of the present paper is to generalize the above results, valid for 
shear-free perfect fluids, to shear-free anisotropic fluids. For this purpose we 
use the reformulation of the Einstein equations in terms of the mass function 
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[IS], [3E]- It is combined with a formula for the same function derived recently 
[12], [E], [IS] and adapted to our goals. In this approach, besides the new 
solutions, the classical results reviewed above follow as particular cases when 
the isotropy condition is imposed. 

In Sec. 2 we give a derivation of the mass decomposition formula in the 
general case of radiating spherically symmetric anisotropic fluid when both 
heat flow and null fluid are present. It consists of a local and a global part, the 
latter being the second Weyl invariant. In Sec. 3 the field equations based on 
the mass function are given for non-radiating spherically symmetric metric. 
The fundamental potential Z is introduced and a new formula for the shear is 
presented. Sec. 4 discusses the general solution of the shear-free anisotropic 
model based on a single second order differential equation. In Sec. 5 the 
particular cases of perfect fluids and charged anisotropic fluids are solved, 
making contact with the existing literature. In Sec. 6 we impose in addition to 
the isotropy condition the requirement for EOS and rederive the Wyman and 
the charged Wyman solution. In Sec. 7 the general conformally flat solution 
for anisotropic shear-free fluids is found and the classic Stephani solution is 
obtained in particular. The general uniform density solution is derived in 
Sec. 8 and the known solution for perfect fluid is reobtained. A mechanism of 
producing models with special linear EOS is outlined in Sec. 9. Sec. 10 studies 
the case when the radial and the tangential pressures are proportional or one 
of them vanishes. Finally, some conclusions are drawn in the last section. 



2 Mass decomposition formula 

Spherically symmetric relativistic fluid models are described by the metric 

ds 2 = e 2 "dt 2 - e 2X dr 2 - Y 2 (d9 2 + sin 2 6d<p 2 ) , (1) 

where is, A and Y are functions of t and r only. The spherical coordinates are 
numbered as x° = t, x l = r, x 2 = 9, x 3 = if. In general relativity one defines 
the Weyl tensor Cijki through the Riemann tensor of the spacetime Rijki as 
follows @9] 

Rijki = Cijki + - (gikRji + gjiRik — gjkRu — guRjk) — - {gik9ji — 9jkgu) R- 

i 

(2) 
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Here R^ is the Ricci tensor and R is its trace, the Ricci scalar. The mass 



function m is given by 
We need the Einstein equations 



m = \yr\ 32 . (3) 



R ij ~ 7;9ijR = %KT ij} (4) 



where is the energy-momentum tensor (EMT) and we have set G = c = 1 . 
For an anisotropic fluid model with heat flow one has [50] 



T ik = (p + p t ) UiU k - p t g ik + (p r - p t ) XiXk + Qi u k + QkUi- (5) 

Here p is the energy density, p T is the radial pressure, p t is the tangential 
pressure, u % is the four-velocity of the fluid, \ % is a un h spacelike vector 
along the radial direction and q l is the heat flux (in the radial direction too). 
They satisfy the relations 

u% = 1, \ l Xi = -1, u l Xi = 0, u% = 0. (6) 

The coordinates are assumed to be comoving, hence 

u^e'^i, X 4 = e-\5i, q* = (7) 

where q = q(r,t). The traces of Eqs (4,5) give respectively 

R = -8vrT, T = p-2pt- p r , (8) 

T being the trace of the EMT. 

Now we are in position to rewrite the basic Eq (2) in terms of Ty instead 
of R^ 

8tt 

Rijki = Cijki — 47T (gjkTu + guTjk — g%kXji — QjiXik) — ^~ {QikQji — 9jk9ii) T. 

(9) 

Let us note that for non-trivial Ui,Xi an d Qii i — or 1- Then for i,j,k,l 
equal to 2 or 3 the previous equation greatly simplifies 

Rijki = C ijM - -£-(p+Pt- Pr) {9ik9ji ~ 9jk9a) ■ (10) 
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This gives for the mass function 



4-7T Y 

m = —Y 3 (p + p t - p r ) + -C 3 232 . (11) 

Note that the tangential pressure enters this formula thanks to the trace T, 
while the heat flow coefficient q does not appear at all. A similar formula 
has been found by inspecting the field equations |47j . where non-comoving 
system of coordinates has been used and the expression in the brackets reads 
Tq +T± — T 2 . In non-comoving coordinates the individual components depend 
also on the radial velocity of the moving observer. One can check, however, 
that in the cited combination this dependence is cancelled and it coincides 
with the bracket term from Eq (11). We have proved this formula without 
resorting to the explicit expressions for the field equations. It is consequence 
of the 'orthogonal' character of the four-vectors present in the EMT (pointing 
in the time and radial direction) with respect to m, defined in terms of tensor 
components along the two spherical angles 9 and (p. 

Formula (11) holds also when the basic anisotropic fluid radiates a null 
fluid with EMT 

If = eltlj, (12) 

where /* is a null vector, l l = u l + \ l • One can show [H] that it adds to the 
basic Tik effective energy density, radial pressure and heat flow, all of them 
equal 

N N N /io\ 

p =p r =q =e. (13) 

This addition does not change Eq (11) or, put in other words, e (like q) does 
not appear in it directly. 

Next we shall use the definition of the second Weyl invariant ^2 in terms 
of the Weyl tensor for the spherically symmetric case [51]. It reads after a 
slight transformation 

2* 2 = S*» = -Y-*C\ 32 . (14) 

5 , 225 l 33 

Hence, one can rewrite formula (11) as 

Att 

m = —Y 3 (p + p t -p r )-Y 3 ^ 2 . (15) 
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In the static case one can set Y — r. For a perfect fluid m contains only 
its energy density. In the general case the anisotropic factor Ap = p t — p r 
is added. Thus, the total mass function is a sum of the local part arising 
directly from matter sources and the Weyl term, responsible for some global 
curvature. This last term may be related also to the electric part of the Weyl 
tensor [IB] 

3 Field equations 

From now on we accept that there is no radiation (q = e = 0). We also 
change the signature of the metric. The usual field equations (4) are rather 
cumbersome except for the (01) component 



The dot above means time derivative, while the prime denotes a radial one. 
In order to make use of the mass decomposition formula (15) we shall work 
in the formalism based on m [15], [46J, namely 



Y' - Yv' - Y'X = 0. 



(16) 



m! = 4npY 2 Y', 
m = -4irp r Y 2 Y, 



(17) 
(18) 



X(p+p r ) = -p 



2Y 



(p + Pt) , 



(19) 



V' (p + Pr) = -P' r + - 

The mass is given by Eq (3) which yields 



Y 
2Y' 



(Pi - Pr) 



(20) 



Y 




(21) 



The expansion of the fluid is 




(22) 



while the components of the shear tensor are proportional to 




(23) 
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Let us introduce next a quantity Z, which is central in the following 



Z = m 2 - —Ap. 
Then the mass decomposition formula (15) may be written as 



(24) 



(25) 



This is an expression for the energy density which, unlike Eq (17), contains 
the mass and not its derivative. Making use of Eq (17) we obtain 



m 



m'Y 
3F 7 



= -Y A Z 



or in a more compact form 



m\> Y 



/ m 



(26) 



(27) 



Next we find another expression for the shear. Replacing m from Eq (15) 
into Eq (18) gives 



3Y(p + Pt )=Yp r -Y(p + p t ) + 



3(Y^ 2 y 
AttY 2 



(28) 



Replacing in Eq (23) A from Eq (19) and using the previous equation we find 



(p + p r ) a 



3Y 



3(Y 3 ^ 2 y 

AttY 2 



- (YAp)' 



(29) 



The formulae in this section hold for any spherically symmetric metric and 
no radiation. 



4 Shear-free anisotropic fluids 

Let us impose on the metric the shear-free condition a = 0. Eq (23) can be 
integrated up to an arbitrary function of r. Usually it is chosen in such a 
way that 
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Then Eq (16) is also integrable and determines v up to an arbitrary function 
of time. Eq (22) shows that it is related to the expansion of the fluid and 
finally 

9 Y 2 

(31) 



-2i< 



@2 y 2 ' 

where G (t) is an arbitrary function pQ . 

In the following we shall use sometimes the notation 



x = r 



(32) 



The metric may be written in terms of Y or L 



ds 2 



9 y s 



Y 



dt 2 + ~—dr 2 + Y 2 dQ 2 



02 y2 
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--^L 2 dt 2 + dr 2 + r 2 dQ 2 



Here dQ denotes the angular part. The same is true for Eq (21) 



9 2 

m = -Y ( 1 + — Y 
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r 2y/2 
~Y2 



Y 



02 

— + 2L X (L - xL x ) 

lo 



(33) 
(34) 



which leads to a formula for the density 



4?r e 2 

— p = — + 2L X (L - xL x ) + Z. 
6 lo 

Eq (27) simplifies considerably because of Eq (34) and becomes 



4 
3 



Z. 



(35) 



(36) 



Let us utilize next the second formula for the shear, namely Eq (29). It 
supplies the equation 



3 (F 3 fr 2 
4vrr 2 



- (YAp)' = 



^2 

L 3 



47T ( Ap 



L — 



(37) 



Combining it with the definition of Z one expresses both Ap and ^ 2 through 
Z 

4vr . (Y 3 Z)- L 4 ( Z 



2Y 2 Y 2L \L 3 



(38) 
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(YZ\ L 2 fZ\ 



2Y 2L \L 

The radial pressure is determined from Eq (18) where the mass is given by Eq 
(34). Together with the above formula for the anisotropy factor and Eq (25) 
for the energy density they give two expressions for the tangential pressure 



Anp t Y 2 Y = - + ^Y 3 ZJ = i (m - AnpY 3 )' . 



(40) 



In the general formalism p and p r are expressed through the mass by the 
similar formulae (17,18) while p t should be found either from Eq (19) or Eq 
(20). For shear- free fluids, however, it is given by the above formula, similar 
in structure to the formula for the radial pressure. 

Now we are in position to outline in two different ways the general solution 
for shear-free anisotropic fluids. First, we choose an arbitrary L(r,t) and 
find the metric, the mass function and Z according to Eqs (33,34,36). Then 
the energy density, the pressures and the second Weyl invariant are given 
correspondingly by Eqs (35,18,40,39). In the second way we take an arbitrary 
Z (r, t) and solve for L from Eq (36), promoting the integration constants to 
arbitrary functions of time. This, however, is rather difficult. For example, 
if Z is a power x n+1 , we get an Emden- Fowler equation which has explicit 
solutions for n = 0, — 1, —2 only [52]. Different additional constraints on Z 
will be discussed in the following, leading to other solutions. 

The shear-free condition (37) can be written in integral form too 

87r f a »,2,> . 9( r ) 87tL3 f a L 



I ApY 2 Ydt + g -^= 8 -^fj Ap^dt + h (x) L 3 (41) 



3Y 3 

where g, h are arbitrary functions of r or i. 



5 Perfect fluids. Charging a fluid 

The arbitrary function in the general shear-free anisotropic fluid solution 
indicates that we are free to impose another condition. This can be done 
in a number of ways. As a first example we impose the isotropy condition 
p r = p t which makes the fluid perfect. Hence, perfect fluid can be viewed 
upon as an anisotropic fluid with linear EOS among the pressures. Then the 
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integral term in Eq (41) disappears and the central equation (36) becomes 



L xx = F[x)L\ F(x) = ^^. (42) 

This equation is well-known [3], [I]. When F = (x — a) n , where a, n are con- 
stants, it becomes an Emden-Fowler equation and possesses five analytical 
solutions for n = 0, -15/7, -20/7, -5/2, -5 [9], [52]. All of them may be ex- 
pressed through the Weierstrass elliptic function. The known solutions of Eq 
(42) PQ, [H], [HI] include some others, like the second Painleve transcendent. 
The time dependence appears in L by turning the integration constants into 
integration functions of time. The general formulae give expressions for the 
other characteristics of the fluid and the freedom in the solution is reduced to 
arbitrary functions of time. Eq (24) shows that for Ap = we have Z = ^ 2 
and hence 

^2 = \xF(x)L\ (43) 

which too is a classical result [I], [16]. We see that the appearance of the 
second Weyl invariant is due to the mass decomposition formula, which in 
its turn stems from the definition of the Weyl tensor, applied to an interior 
fluid solution. 

Another way to constrain the form of Ap is to charge the fluid. Spherical 
symmetry allows the appearance of only an electric field E in the radial 
direction. The energy-momentum tensor of this field can be described [H] 
as addition of effective energy density and pressures to the fluid with 

p E = pf = ~pf = f- (44) 
The Maxwell equations give 

Ap E = ^ = / (x) L\ 4vrr = Ee~\ (45) 
4tt 

where r is the charge function of the fluid and / (x) is an arbitrary func- 
tion. Hence, Eq (41) acquires a charge term in addition to other sources of 
anisotropy 

Z = h(x) L 3 + ^-f (x) L 4 + |^L 3 f Ap^dt. (46) 
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In the case of charged perfect fluid Eq (36) becomes 

xx = J~x L + "T" L ' (4?) 

which is also well-known and studied in the past [TT] . [12]. The Weyl invariant 
reads 

^ 2 = hL 3 + 4vr/L 4 (48) 
and we have shown the reason for its close resemblance to Z. When 

h = x ni+ \ f = x n2+1 (49) 

Eq (47) has five different analytical solutions for (ni, n 2 ) = (0, 0) , (—5/2, —3) , 
(-14/5, -18/5) , (-11/5, -12/5) , (-5, -5) [52]. The density and the pres- 
sures are given by the general formulae plus the additions from Eq (44). In 
total, charged perfect fluid is a particular case of the anisotropic fluid model. 

One can find new neutral anisotropic solutions by prescribing the form of 
the anisotropy factor. For example 

Ap = f(x)L a , L xx = ^I? + -^—L a -\ (50) 

4x (a — 3) x 

where a is some constant, not equal to 3. 



6 Solutions with equation of state 

There is a case when the integration of Eq (47) is straightforward. Let the 
x-functions on the r.h.s become constants a,/3. Multiplication with L x gives 
a first integral 

Ll = e + l*L 3 + \PL* (51) 
where £ is another constant. Integrating once more one finds 



/ 



x + t = u, (52) 



| L 4 + + £ 



where an arbitrary function of time has been set to t. This is a particular 
form of the integral 

u= f (s) -1 ^ 2 ds, f (s) = aos A + Aais 3 + 6a 2 s 2 + 4a 3 s + a 4 (53) 
Jo 
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and other examples of it will appear later. It may be inverted in terms of 
the Weierstrass function P (u, #2, #3) [S3], namely 

2(P-f) 2 -ia a 4 ' 

(72 = flo a 4 — 4axa 3 + 3a 2 ,, (55) 
#3 = a ct2O4 + 2aia 2 a3 — — ao a 3 — a 2 ^. (56) 
Here ai are constants or functions of time. In the present case 

92 = g 3 = -^« 2 e (57) 

Eq (35) shows that while L = L{u) the density in general is p (t, r) because 
of terms with x. However, Eqs (44,45) lead to their cancellation and tuning 
the expansion yields 

%ixp (u) = 3 (C - 4£u) + 12LL U , 6 2 = 9C - 36£t, (58) 

where C is a constant. Eq (19) holds for the effective energy and pressures 
of the charged perfect fluid. It gives the following equation for the sum of 
the genuine energy density and isotropic pressure 

47-4 /r o \ 
Stt (p + p) = — l-a+lpL). (59) 

Hence, p = p(u) and p = p{p). The fluid possesses an EOS and this is the 
charged Wyman solution [23] . [24] . [25] . [26J. Combining the two previous 
equations we obtain 

L L 5 

8tt (p + 6p) = -15 (C - A^u) - 60£— + 6/3—. (60) 

When the fluid is neutral (^ = 0, (3 = 0) the EOS becomes linear and the 
Wyman solution emerges [3], [22] . 

In principle, the general solution of the shear-free anisotropic fluid de- 
pends on an arbitrary function L(t, r). Going to a perfect or charged perfect 
fluid means imposing a differential equation on L so that the freedom reduces 
to functions of t or r. Imposing further an EOS makes the system overdeter- 
mined and solutions (if any) may exist only in some exceptional cases, like 
Eq (60), which holds for constant £, a, (3. 
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Is it possible to generalize this EOS mechanism to anisotropy which is 
not the result of charging the fluid? For this purpose we take 

Z=^xLF(L) L , L XX = F(L) L , L 2 X = 2F (L) + g (t) , (61) 

where F (L) is an arbitrary function of L. The main equation still may be 
integrated and L = L (u) if the integration function of time is constant go- 
Eq (35) gives 

e 2 

Aixp = — + 6LL U + 2x (2LF L - 6F - 3g ) . (62) 

The elimination of the x-term requires an equation for F leading to the 
Wyman solution. Thus, it is not possible to generalize it to anisotropic 
fluids except for a charged perfect fluid. 



7 Conformally flat solutions 

These solutions require * 2 — 0. The definition of Z and the vanishing shear 
constraint (37) then yield 

Ap = f(x)L, Z = ~f(x)L, L„ = -^, (63) 

O Ob 

f (x) being an arbitrary function. The differential equation is easily inte- 
grated so that 

L — G(x) + h(t)x + q(t), f (x) = --xG (x) xx (64) 

71 

and L depends on three arbitrary functions. The expressions for m,p,p r ,p t 
are rather cumbersome and are not given. Let us note a general equation 
which follows from Eqs (18,40) 

L . . 

2p t + p r + 3p = j p. (65) 

Things simplify considerably for a conformally flat perfect fluid. Then 
one can put G = and Eqs (15, 35, 64) give 

9 it) 2 4irr 3 o 
Anp = ^L + Qh(t)q(t), ™ = ^j/, (66) 
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so that p is a function of time only. Eq (65) then gives for p 



V=- P +'f L - (67) 

This is exactly the Stephani solution PQ, [36]. Eq (29) shows that when 
^>2 = = Ap the shear vanishes identically. Thus Stephani's solution de- 
scribes all conformally flat perfect fluids, while the solution (64) describes all 
conformally flat anisotropic fluids with vanishing shear. Comparison between 
Eq (45) and Eq (63) shows that the anisotropy can not be caused by charging 
the fluid unless the solution is static, hence, evolving charged perfect fluid 
can not be conformally flat in addition. This result was found in the past by 
other methods 



8 Uniform density 

In the previous section we have seen that some conformally flat solutions 
possess uniform density p (t). One can find all solutions of this kind. Up to 
now we have worked in the L-formalism, taking L as a basic function. From 
now on we shall work in the ^-formalism. For uniform energy density Eq 
(17) may be integrated 

Air 

m = —pY 3 - f (t) (68) 

and Eqs (27, 38) give 

Z = M, ^Ap=— L. (69) 
F 3 3 y 2Y 2 Y 

The combination of Eqs (34, 68) provides a first order differential equation 
for Y 

If = AY^ + Y 2 + 2fY, (70) 

a V) = y~T p ' 2: = lnr ' ( 71 ) 

which can be integrated 

I , dX = = s, X = l s = \nrB(t). (72) 
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Here another arbitrary function of time appears. One can invert the func- 
tional dependence with the help of formula (54) and express Y in terms of 
the Weierstrass function with 

1 1 f 2 A 

92 = T2 , a. = - - — ( 73 ) 

The pressures are obtained from Eqs (65, 69). 

In the case of perfect fluid the function / (t) becomes a constant and the 
above solution still applies. When this constant vanishes, Z = = ty 2 and 
the solution is conformally flat. The integral above is expressed in elementary 
functions. We find 

rX = L = e- x = -- + -r 2 , (74) 
B 4 ' v ' 

which coincides with Eq (64) when G — 0. These results are the same as the 
classical results on perfect fluids with uniform density [I], [27]. They can be 
rewritten in the L-formalism. There L satisfies Eq (42) with F (x) ~ x~ 5 ^ 2 
which follows from Eq (69). 

In the general anisotropic case one has Y = Y (s) when / and A are 
constants. Then Ap = and the fluid is in fact perfect. Further, m = m(s) 
when is constant. This leads to the constancy of p and p. Then Eqs (18, 
68) show that p+p = 0. Once again the fluid has an EOS when the arbitrary 
time functions become constants. 



9 Solutions with p + p r = 

We may generalize the uniform density case with the help of an arbitrary 
function F (Y). Mimicking Eq (69), let us take 

Z = f(t)F(Y) Y Y (75) 

In particular, uniform density follows from F = — 1/3Y 3 . Eq (27) gives the 
mass function 

Tfl 

— 3 =3f (t) F(Y) + g (t) (76) 
and then p is found from Eq (25) 

47T 

—p = f(3F + YF Y ) + g. (77) 
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Proceeding like in the previous section we obtain 



dX 



A x (*) = — - 2g, 



(78) 



yJ-efF(i)+X^ + A 1 



which determines Y (s,t). We also find 



4tt (p + p r ) = - (3fF + g) -^-Ap =fF Y ^ + 4Z, (79) 



which determine the pressures. When / (t) ,g (t) are constants we get a 
linear EOS p + p r = 0. The density and pressures are not constant like in 
the previous section. 

The appearance of EOS can be seen in another way too. We have m = 
m(Y) from Eq (76). Eqs (17,18) indicate that in this case the difference 
between time and space derivatives effectively disappears and p + p r = 0. 
The discussion in the present section is analogous to that in Sec. 6 which 
lead to the Wyman solution. However, now we work in the ^-formalism and 
the EOS differs from Eq (60), reached in the L-formalism. 

10 Solutions with p t = or p r = 

Solutions with vanishing tangential or radial pressure have been discussed in 
the past for time-dependent anisotropic fluids with shear [51], [55], [56]. We 
shall investigate the more general case of linear EOS between the pressures 



Here 7 is some constant measuring the degree of anisotropy between the two 
extreme cases 7 = (p t = 0) and 7 = 1 (pt — p r , isotropy condition). We 
have shown that for shear-free fluids pt satisfies Eq (40), similar in structure 
to Eq (18) for p r . This allows to write Eq (80) as a time derivative and to 
integrate it, obtaining 



g (r) being arbitrary. In addition, Z satisfies Eq (26) for any spherically 
symmetric metric. This leads to 



Pt = 7Pr- 



(80) 



^Y 3 Z + (l- 1 )m = g(r), 



(81) 




(82) 
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Suppose that g (r) = g , a constant. Then another integration gives 

m = a + k(t)Y 1+2 \ a = (83) 

and (t) is another arbitrary function. Eqs (81, 17) read 

Z = -^ + 2 -^kY^\ (84) 
4vrp= (l + 2'y)kY 2 ^- 1 \ (85) 

Eqs (38, 80) result in 

. y2 7 -i 

47rp r = -Aixp - fc , (86) 

which allows to find \&2 from the definition of Z 

n 1 1 y 2 7-i 

*2 = + 3 (7 " 1) (1 " 27) ^ 2(7 - 1} - 3 (7 - 1) fc— ■ (87) 

The combination of Eqs (34, 83) provides a first-order equation for Y or 
X — 1/Y which can be integrated 

dX , s 

s. (88) 



y/-2aX 3 + X 2 - 2kX 2 ( l ~^ + 2 /9 



The integral may be inverted with the Weierstrass function when the poly- 
nomial Q under the root is up to the fourth degree. This happens when 
7 = 0,±1/2,±1. 

In the case 7 = (pt — 0) the polynomial and the other characteristics 



of the fluid read 



e 2 

Q = -2aX 3 + (1 - 2k) X 2 + — , (89) 

9 

m = a + kY, Z = ~-^ (90) 

4vrp = ^, 4vr(p + p r ) = -i, (91) 
T ci k k . . 

^2 = -TT^ + 7TT77 + 92 

F 3 3F 2 3FF 
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When k is constant p + p r = 0,m = m (Y), Z = Z (Y) and the solution 
belongs also to the class discussed in the previous section. 
In the case 7 = 1 (perfect fluid) one finds 



Q 



-2aX 3 + X 2 + — -2k, 



(93) 



a + kY 3 , 



Z 



a 



(94) 



m 



y3 ' 



Anp = 3k, 



4tt (p + p r ) 



kY 

T' 




(95) 



Now p = p{t) and this is in fact a uniform density solution given by Eq (72) 
when / it) = —a. When a = the solution is conformally flat as can readily 
be seen from Eq (95). 

Let us discuss finally the case p r = 0. Eq (18) gives m = m(z). Substi- 
tuting this into Eq (34) yields 



This equation may be integrated easily when m is constant and a formula 
close to Eq (72) or to Eq (88) with Q from (89) or (93) is obtained. However, 
this is unphysical case because the energy density vanishes. Therefore Eq (96) 
should be solved by other methods. 

11 Conclusions 

Recently we have stressed the importance of anisotropic fluid models for 
astrophysics [H]. With the present paper we start their systematic study 
from the most general case. A decomposition of the mass function into a 
local and global part is redireved and related to the second Weyl invariant. 
It holds for any radiating anisotropic fluid which is spherically symmetric. 
Then, for simplicity, we turn off the radiation and combine this mass formula 
with the reformulation of the Einstein equations in terms of the mass function 
(given by a classic formula) and its derivatives. This results in the important 
relation (27) for the master potential Z and another expression for the shear, 
Eq (29). Here the usual metric formula for the shear (23) is elevated to global 
and source quantities such as \&2 and the anisotropy factor. It allows to push 



X 2 = -2m (z) X 3 + X 2 + 



9 



(96) 
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further the study of shear- free fluids, providing an important constraint. Now 
everything can be expressed in terms of L and Eq (27) transforms into the 
second order Eq (36). The general solution of this equation is based on one 
arbitrary function which may be Y, L, Z, Ap or ^2- The treatment of the 
tangential pressure becomes on an equal footing with the radial pressure due 
to formula (40). The acceptance of maximum freedom in the system allows 
to obtain an anisotropic shear-free core of simple relations. Only after that 
we start to impose one by one different constraints aiming the system at 
particular cases. 

Thus perfect fluid is this core with imposed isotropy condition (Ap = 0), 
charged fluid is a neutral fluid with a special kind of anisotropy, conformal 
flatness means the constraint ^2 = 0, uniform density constrains the energy 
density to a function of time, anisotropy may be studied in more detail by a 
linear equation among the pressures, which includes the cases where one of 
them vanishes. In this way the relations among the particular cases become 
much clearer due to the relations of all of them to the core. Time and 
again second order equations like Eq (50) appear, which sometimes may be 
integrated to first order, like Eq (61) and even to implicit integral formulas 
for the solution, like Eqs (52, 72, 78, 88). These are different examples of 
the integral (53) which explains the persistent appearance of the Weierstrass 
elliptic function. In many cases, however, it reduces to elementary functions. 
This is true also for the ^-formalism. 

Imposing one constraint fixes the arbitrary potential of the general model 
but up to functions of time or radius. Imposing a second constraint makes 
the system overdetermined, yet special solutions still exist when the arbitrary 
functions become constants. Then the integral formulas reduce the depen- 
dence on t and r to a dependence on a single variable r 2 + 1 or rB (t). The 
result is appearance of equation of state, linear in some cases, as seen from 
Eqs (60, 79). 

In total, we find the general charged anisotropic solution, all conformally 
flat solutions, all uniform density solutions, a large class of solutions with pro- 
portional or vanishing pressures and solutions with vanishing sum of density 
and radial pressure. 

Finally, let us discuss the static case. Unfortunately, Eqs (16-21) are 
incomplete then because (16,18,19) are trivial. The shear trivially vanishes 
and formula (31) holds no more. The metric function v is not determined by 
Y and becomes another free potential. One can use it and Ap to express all 
static spherically symmetric anisotropic solutions [57J. 
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